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O ■ Abstract 

| — 

consequence amplification of gauge fields. The resulting spectrum of the field strength is ap- 
proximately Bt oc i" 1 , where £ is the relevant coherence scale. One realisation of our scenario 
is scalar electrodynamics with a scalar whose mass is large enough to evade observational 
^ ■ constraints - the obvious candidates being supersymmetric partners of the standard-model 
fermions. Our mechanism also leads naturally to amplification of the standard-model Z- 
boson field due to its coupling to the electroweak Higgs field. At preheating, the spectrum 
of the Z field is transferred to the hypercharge field, which remains frozen in the plasma and 
is converted into a magnetic field at the electroweak phase transition. With a reasonable 
model of field evolution one obtains a magnetic field strength of the order of 10~ 29 Gauss 
on a scale of 100 pc, the size of the largest turbulent eddy in a virialised galaxy. Resonant 
amplification in preheating can lead to primordial fields as large as 10 -24 Gauss, consistent 
with the seed field required for the galactic dynamo mechanism. 
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1 Introduction 



It is well known that during inflation the inflaton, being a light scalar field, couples gravi- 
tationally and grows on superhorizon scales. Likewise, (scalar) cosmological perturbations 
grow resulting in a scale-invariant spectrum on superhorizon scales, providing one of the 
most important predictions of inflation: seeds for formation of large-scale structures in the 
Universe fl], § . The spectrum of gravitational waves generated during inflation has also been 
studied in detail and is regarded as perhaps the only signature of inflation that may be 
"observed" directly. It is hence surprising that the evolution of vector fields has not been 
considered with sufficient care, especially in light of the exciting possibility that primordial 
gauge fields from inflation could produce the magnetic fields observed today in galaxies Q 
and even affect the cosmic microwave background radiation 

It is usually assumed that gauge fields do not grow during inflation. The reason is very 
simple: gauge fields do not couple gravitationally to a conformally flat space-time. The 
metric of conformally flat space-times can in general be written as g^ v = a 2 (r, x)^^, where 
Vtiv = diag[l, —1,-1,-1] is the Minkowski metric. The space-times of standard cosmology 
are all conformally flat. Indeed, in inflationary de Sitter space-time the scale factor reads 
a(r) = —1/ (Hit), where Hi is the Hubble parameter in inflation; in the radiation era a oc r, 
while in the matter era a oc r 2 . This implies that, in order to get any amplification of gauge 
fields during inflation, conformal invariance must be broken. 

Various proposals have been put forward in which conformal invariance of gauge fields is 
broken. A particularly nice one is due to Dolgov || , who showed that the conformal anomaly 
of the gauge-field stress-energy tensor, expressed by the triangle diagram, leads to particle 
production for gauge fields even in a conformally flat space-time. Turner and Widrow || 
consider several ways of breaking conformal invariance: (a) gravitational coupling of the 
photon, (b) anomalous coupling of the photon to the axion, and (c) coupling of the photon 
to a charged, massless scalar field. In the former two cases, conformal anomaly is broken 
by a non-standard coupling of the photon field. The case (c) has recently been explored by 
Calzetta et al. and by Kandus et al. 0. These authors argue that, as a consequence of scalar 
charge separation during inflation, charged domains form which proceed to source electric 
currents in the radiation era, leading to magnetic-field production. Their claim, that the 
magnetic field thus created has sufficient strength to seed the galactic dynamo mechanism 
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||, has recently been contested by Giovannini and Shaposhnikov 0. 

In this Letter we show that there is a natural way of breaking conformal invariance of 
gauge fields which has not been considered in the literature. Namely, the backreaction of a 
charged scalar field gives the gauge field an effective mass, which breaks conformal invariance. 

2 Scalar electrodynamics 

We begin by considering scalar electrodynamics, which nicely illustrates our mechanism of 
gauge-field amplification from inflation. The idea can be easily extended to include non- 
Abelian gauge fields and multi-component scalar fields. The Lagrangian is 

C^d = -\g m g up F^F vp + g^(D^)W u - mj0t _ A (0t0) 2 , (1) 

where = — ieA^ is the covariant derivative, F pv = d^A v — d u A p is the gauge field 
strength, g^ u = a~ 2 r]' J ' u is a conformally flat metric, and and A^, are the mass and 
self-coupling of the scalar field 0, respectively. The experimental constraint on the mass, 

3> Mew — 100 GeV, is easily satisfied. The obvious candidates for are supersymmetric 
partners of the standard-model leptons and quarks. In addition, we require that <C Hj, 
so that the charged scalar field may grow during inflation. 

The relevant equation of motion for the photon field is obtained from Eq. ([!]) as the 
Lagrange equation d v (5g vA)1 ( ^-DgC^v)) ~ S Afl ( ^-A^ed) = 0, where D g = det^] = 
—a 8 . After some algebra one finds that the mode equation for the transverse component of 
the photon field in the Hartree approximation and in unitary gauge can be recast as 

(^ + £ 2 + e 2 aV>)^ = 0, (2) 

where p 2 = 20^0, k is the comoving momentum, e electric charge (e 2 /47r = a ~ 1/137), 
and r the conformal time, which is related to the "physical" time t as dt = a(r)dT. The 
average (•) should be computed in the Hartree approximation by subtracting the vacuum 
contribution, so that in the vacuum (p 2 ) — > and Eq. (|2]) reduces to that of a harmonic 
oscillator. The solutions are then travelling harmonic waves, 

.e> = (^)v-, ( 3 ) 
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where k = \k\ and we used the (Wronskian) normalisation of the mode functions, W[^44+j, A^. 
= AiriV, appropriate for periodic boundary conditions in a box of volume V. The modulus of 
the mode functions in Eq. (^) does not depend on time, which is a consequence of conformal 
invariance being restored when (p 2 ) — > (p 2 ) vac = 0. 

For simplicity we shall consider only homogeneous space-times such that the scale factor 
a in Eq. (§) may be taken to be 

_ f -I /{Hit) for r < -H{ 1 (inflation), 
1 H\t for r > Hi (radiation). 

This corresponds to the following smooth matching at the inflation-radiation transition: 

a(-H{ v ) = a(Hr l ) and da(-Hf l )/dr = da(H{ 1 ) / dr . 

2.1 Breakdown of conformal invariance 

The backreaction term in Eq. (|2|) is generated during inflation because of the gravitational 
coupling of the scalar field. In the Hartree approximation, the backreaction is described by 
elastic scattering processes only, which should be a good approximation in inflation. We 
also assume that damping is negligible in inflation, since most of the amplification is on 
superhorizon scales beyond the reach of dissipative processes. 

The Hartree term (p 2 ) for a light, charged scalar field can be estimated as follows. From 
the scalar equation of motion, one infers that during inflation the scalar field grows until 
its effective mass essentially reaches the Hubble parameter. More precisely, we have m\ + 
3A (p 2 ) < 2H\, so that e 2 (p 2 )/H? < 2e 2 /(3X (j> ) ~ 1/(16A^). When the limit is saturated, 
the growth of 4> modes becomes exponentially suppressed. We assume that this occurs well 
before the end of inflation. The solutions of Eq. (0) may then be conveniently written in 
terms of Hankel functions as follows: 

Af =Tc(-TV)*H®(-kr), t = l,2, ^ 2 = I-^1 (inflation). (5) 

For the mode-function normalisation we impose the Wronskian condition W ^4.^ , = 
AiriV so that, as r — >• — oo, the mode functions reduce to the vacuum mode functions 



I (1,2) 
k 

get "squeezed" on superhorizon scales as 



discussed above: A// — > -4^. At later stages of inflation, however, the mode functions 



Since Re [u] < 1/2, the amplitudes A~ actually decrease, while the momenta d T A^ grow 
(see Fig. |I|). This evolution preserves the uncertainty relation for the photon field operators, 
which can be equivalently expressed by the (conserved) Wronskian. 

After inflation the squeezed photon state evolves according to complicated physics asso- 
ciated with the processes of preheating and thermalisation of the Universe. The details of 
the evolution are model dependent and deserve further study [pj, especially since these pro- 



cesses have not been fully considered in the context of gauge fields (see, however, Ref. [I I]]). 
In the following section we shall make a crude approximation and assume that the scalar field 
decays instantaneously at the end of inflation, leading to sudden restoration of conformal 
symmetry. The effects of preheating and thermalisation are considered in section |, where 
we show that, under certain conditions, preheating results in additional amplification of the 
superhorizon modes. 

3 Photon field amplification from inflation 

Here we assume that soon after inflation the <fi field decays nonadiabatically such that (p 2 ) 
quickly approaches the vacuum value (p 2 ) vac = 0, rendering Eq. (E3) conformally invariant. 
The solutions in the radiation era are then linear combinations of the travelling waves A^ 
of Eq. ([|). Particle production occurs since the effective mass parameter undergoes a non- 
adiabatic change at the inflation-radiation transition. This change can be nonadiabatic even 
when the squared mass e 2 (p 2 ) changes adiabatically. Indeed, at the beginning of radiation 
era the scale factor evolves nonadiabatically, i.e. \(d T a 2 )/a 2 \ 2 = AH 2 ^> e 2 (p 2 ), so that the 
results concerning the photon field amplification in the simple case under consideration are 
quite generic. 

Since the vacuum of inflation does not correspond to the late-time vacuum of the radiation 
era, the field content of Eq. ([5]) may be interpreted from the point of view of a late-time 

,A {1) 

combination of travelling waves in the radiation era, 



observer by a smooth matching of the mode functions A~ at the end of inflation to a linear 



drAfi-Hr 1 ) = a^AfiH^ + ^drA^Hr 1 ), (7) 
and similarly for A^ = A^* . Here and are the Bogoliubov coefficients that relate 
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Figure 1: Evolution of (the imaginary part of) the gauge-field amplitude Im [At] (red solid curve) 
and its momentum Im [<9 T .Ag] (blue dotted curve) for k/H\ = 0.05 near the inflation-radiation 
transition, which occurs at r = r, = ±H7 . As a consequence of broken conformal invariance, 
toward the end of inflation the amplitude decreases and the derivative increases, leading to an 
enhanced amplitude in the radiation era. The normalisation is chosen so that both functions would 
have amplitude equal to unity in the conformally evolving case, and we use v = 0.2 for the purpose 
of illustration. 

the creation and annihilation operators in the radiation era to those of (de Sitter) inflation!] 
With A^ = A^ the matching problem is easily solved. On superhorizon scales we have 

so that at late time in radiation the linearly independent mode functions A^ and At*, where 

Aj: = a^A^ + P*t^_t > are proportional to k~ x ~ v . When v « 1/2 the superhorizon modes 

«4.£ oc fc -3 / 2 make up an almost scale-invariant spectrum, similar to that of a massless scalar 

field. The amplification by a factor (3^ oc {kj H\)~ x l 2 ~ v compared to the Minkowski vacuum 

spectrum (see Fig. [I]) is a consequence of the large value of the momentum d T A*jj at the end 

of inflation, as can be inferred from Eq. (^|). This is the maximum amplification that can be 

attained at the inflation-radiation transition, since sudden transition is the one with maximal 

nonadiabaticity. However, it is not hard to show that a softer matching, for example onto 

The convention used here is consistent with the Bogoliubov coefficients being defined through the follow- 
ing relations: — a^aj: + — a^b^ + /%5' -., where barred operators are for inflation and unbarred 

for radiation era. 
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(p 2 ) ~ H 2 = (Hit 2 ) 2 , leads to almost as strong an amplification on superhorizon scales (for 
v ~ 1/2), namely At oc k~^~ 2v . An obvious candidate for such a scenario is "gravitational" 



preheating [0. We emphasise that, in a transition between asymptotic "in" and "out" 
states with the property {p 2 )- in = (p 2 )out = 0, the expectation value (p 2 ) will typically be 
non-zero during any period when the quantum operator p 2 evolves non-adiabatically, leading 
to amplification of gauge fields. 

4 Photon field amplification in preheating 

Here we present a simplified model of preheating in which the spectrum of photons from 
inflation is further amplified on superhorizon scales, so long as inelastic scattering and dis- 
sipation can be neglected. We assume that, at late stages when thermalisation takes place, 
the (j) field decays nonadiabatically, imprinting the photon-field spectrum onto the plasma 
as described in section |[ How realistic this assumption is will be the subject of a separate 
study [0. 

Resonant amplification of the photon field may result when the scalar field <fi couples to the 
inflaton and grows resonantly such that the photon mass term e 2 ((p 2 ) (<p 2 = 2a 2 </>V) acquires 
an oscillatory component, leading to (secondary) resonant amplification of the photon field 
A. This resonant amplification resembles stochastic resonance [|13| in which all superhorizon 



modes are equally amplified. We shall now make a crude estimate of the amplification factor. 

We assume massless chaotic inflation with a quartic interaction term A s s 4 /4, where X s ~ 
10~ 13 as specified by the COBE satellite microwave background radiation measurements. 
Furthermore, the inflaton s couples to <fi with a term hs 2 (ft(f) such that hs 2 < 2H 2 (s) during 
inflation, setting an upper limit on h. It is now easy to show that the oscillating inflaton 
s decays into <ft through a narrow resonance with the quality factor « hs^/Auj 2 < 1/200, 
where so ~ 0.3Mp is the inflaton amplitude at the end of inflation and uj\ = c n \/X^ Sq 
(c n ~ 0.847) is the inflaton frequency. Since at very early stages the photon-field quality 
factor qA = e 2 '(f 2 ) osc / '4cof q</>, where (<p 2 ) OS c is the amplitude of the oscillatory component 
of (<p 2 ), energy density p is much more efficiently transferred from s to <fi than it is from <p to 
A. As time passes, {<p 2 ) OS c grows, leading to more efficient resonant production of the photon 
field A. Quite generically, the photon-field superhorizon modes begin to grow when qA ~ 1 
is attained [|13|]. At that time n, Pa(ti) < P4>( r i) ~ 12A^u;^ / e 4 . The modes stop growing 
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when p<£ ~ pa ~ p s — w i/(4A a c^). On the other hand, a resonating vector mode grows as 



Aj: oc e^ AUJlT and hence p^ oc e 2fJ,AlJJlT , implying the following estimate of the amplification 
factor x f° r superhorizon modes during preheating 

X ~ f — 7~T ] 2 ~ 7T~?== ~ 10 5 . (9) 



We assume that the photon field at this point freezes in due to large conductivity. The 
purpose of this simple analysis was to show that superhorizon modes are amplified during 
preheating, a rough estimate of the amplification factor being Eq. (|SJ). To obtain a more 
accurate value, the problem should be reanalysed using numerical techniques. We have also 
studied the evolution of the spectrum by matching the inflation mode functions to those 
in preheating, A^ h = (2ttV/uji) 1 ^ 2 e^ ATt " >u)lT . At the end of preheating, the spectrum is 
v4.p re h oc k~ u , as in Eq. (H), with an additional amplification factor given by Eq. @. A 
final matching of the preheating photon modes onto the massless photon modes in radiation 
results in the spectrum ^4 ra d oc k~ x ~ v computed in section §, again with the additional 
amplification factor (|9|). 

5 Magnetic field production in the standard model 

In this section we show that the amplification mechanism, which we illustrated with the 
example of scalar electrodynamics, is operative for the standard-model Z field just as it is 
for any gauge field that couples to a light scalar field. Namely, in the Z-field mode equation 
{of- Eq. ®) 

conformal invariance is broken by the Hartree backreaction term of the standard-model 
Higgs field $. Here v = 246 GeV and p 2 = 2$ t $. The restriction 3\ H (p 2 } < 2Hf im- 
plies that, in order to obtain substantial amplification of superhorizon modes {y >, 1/3), 
we require \h ^ 0.66. Equivalently, the Higgs-boson mass must satisfy the lower bound 
ran = v \/2\ H > 280 GeV. The spectrum reads 



,2 



dl + e + ^a\p 2 ) Z £ = (10) 
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Figure 2: Projection of the primordial spectrum of the standard-model Z field first onto the 
hypercharge field Y (at the inflation-radiation transition) and then onto the photon field A (at the 
electroweak phase transition). 

where zf ] = 2^ = (2nV/k) 1/2 e Tlkr (cf. Eq. ©), and a % and /% are given in Eq. (§). 
The non-Abelian component of the Z field, i.e. develops a magnetic mass and becomes 

screened as the Universe thermalises. Only its Abelian component, the hypercharge field Y = 

— sin 9yjZ, survives and freezes into the plasma. Here #w is the Weinberg angle (sin^w ~ 
0.23). After the electroweak transition it is the photon field which remains unscreened 
with amplitude A = cosOyjY. Note that, although the Z field and the photon field are 
orthogonal, the resulting photon spectrum is identical to that of the original Z field, apart 
from an amplitude suppression factor sin^wcos^w ~ 0.42. 

The situation is analogous to that of light polarisers (see Fig. |2|). With two orthogonal 
light polarisers no light passes. When a third polariser is inserted at an angle 9 (with respect 
to the second polariser) however, some of the photons do pass. The photon amplitude is 
reduced by sin 9 cos 9, just as in the Z-field case. The main advantage of the amplification 
mechanism presented here is its naturalness. Indeed, no fields are required except those of 
the standard model and inflation. 

One should also consider the possibility of additional enhancement of the field Y = 

— sin^w-^ on superhorizon scales during preheating. For this purpose, we envision a hybrid- 
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inflation model in which inflation is driven by the vacuum energy of a GUT Higgs field. 
After inflation ends, this field begins to oscillate, leading to resonant production of some 
superheavy gauge fields Xj. As a consequence the hypercharge field Y, which couples to 
the Xi, grows as well, and freezes in when the Universe thermalises. The details of this 
amplification mechanism are model-specific and will be presented elsewhere [|ItJ . Finally, we 
emphasise that the superheavy gauge fields grow resonantly during preheating and suddenly 
decay after preheating when they become supermassive. These are required properties of 
the scalar field in the scalar-electrodynamics model considered in sections || and |[ 

6 Magnetic field spectrum 

We now consider the possibility that primordial gauge fields from inflation can provide the 
necessary seed for the galactic dynamo mechanism || and thereby explain the observed 
presence of micro-Gauss-strength magnetic fields in a large number of spiral galaxies ||. 

The dynamo mechanism amplifies a weak, coherent seed field using the differential rota- 
tion of a galaxy in conjunction with the turbulent motion of ionised gas. The length scale 
relevant for the operation of the galactic dynamo is the size of the largest turbulent eddy, 
~ 100 pc. On this scale, there is a minimum strength of the magnetic seed field for which 
the dynamo can operate. Estimates of this minimal strength lie in the range 10 -23 -10~ 19 G 
for a universe with critical matter density and zero cosmological constant. The lower bound 
can be relaxed |14| to about 10 -30 G for a flat, low-density universe with a dark-energy 
component (e.g. a cosmological constant or quintessence), which appears to be favoured by 
recent results from supernova observations and balloon experiments |15L jl6f. 



The dynamo scale (100 pc) corresponds to a comoving scale today of £ c ~ 10 kpc ~ 
1.56 x 10 36 GeV -1 before gravitational collapse and galaxy formation fH |. Because the 



concentration of matter into a galaxy brings about an amplification of magnetic fields by 
a factor (p ga i/po) 2 ^ 3 ~ 5 x 10 3 , the bounds that should be imposed on the scale £ c are 
-Bseed ~ 2 x 10~ 27 G for a universe with critical matter density and -B S ccd ^2x 10 -34 G for 
a flat, dark-energy dominated, low-density universe. 

We consider first the spectrum resulting from scalar electrodynamics (section |3|) and the 
Z-field case (section ^) and neglect for the time being further amplification which may re- 
sult from parametric resonance during preheating. At the end of inflation, the scale factor 



9 



corresponds to an equivalent temperature T given by Hi = ngl /2 (T)T 2 /(V90M P ), where 
g*{T) ~ 10 2 is the number of relativistic degrees of freedom. Assuming first that the primor- 
dial magnetic field is frozen in the plasma from the time of its creation, the relevant length 
scale for the magnetic field at the end of inflation is £ = £ C T Q /T, where T ~ 2.73 K is the 
temperature today. 

For k = 2tt/£ <C Hj, the Bogoliubov coefficients and f}^ given by Eq. (R) are maximal 
for v = v 1/2 — 1/ \n{2H\/k) . For this optimal value of v one obtains the mode function 



M~?W- 2 [ Wl ) ■ (12) 

The magnetic field correlated on a particular scale £ is defined by 

Bj = {Bi{i, x)B { (£, x)} - (B t (£, x)B { (£, x )) vac , (13) 
5 ^> f ) = A/ , *VBM> ( 14 ) 

A-k£ 6 J\y-x\<£ 

where the average (•) is taken over Fock space as well as the position x. For £H\ ^> 1 one 
finds 

Be = 3x 2 p - 2 ^(i - v)H?(£Hy-l . (15) 

7T 2 

The magnetic field then is assumed to be frozen in the plasma such that today Be c = 
Bi(T /T) 2 . Since roughly B^ c oc T, both B^ c and Bg are larger for higher inflation scale 
Hi oc T 2 . Taking Hj ~ 10 13 GeV (corresponding to T ~ 10 15 GeV) we find B t ~ 10 22 G at 
the end of inflation, and B^ c ~ 10 -34 G today at a comoving scale of 10 kpc (see Fig. |3[). 

In section [| we showed that this field strength can be further amplified by a factor of about 
10 5 through parametric resonance with an oscillating scalar field, leading to Bg c ~ 10~ 29 G. A 
supplementary increase in field strength is obtained if we assume that the magnetic field does 
not freeze into the plasma upon creation, but rather that its correlation length grows quicker 



than the scale factor, as is the case for helical turbulence |17| . Such a causal mechanism can 
only operate on a given comoving scale after this scale has reentered the horizon. One can 
show that the growth of correlations due to turbulent evolution leads to an additional am- 
plification of about (£ eq / £ c ) 2 ^ u ^ 3 , where £ eq ~ 50 Mpc denotes the equal matter-radiation 
horizon today. For £ c = 10 kpc the amplification is about 20. For other types of turbulence 
the amplification may be somewhat smaller. 
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(critical universe) 



_ 2x10 Gauss 




-34 

2x10 Gauss 



/ H =H,T/r =10 pc 



Figure 3: Magnetic-field spectra and relevant seed-field bounds. In green (dash- dot- dot- dot) we 
show the vacuum spectrum Bg oc i~ 2 obtained from preheating, assuming an amplification factor 
of 10 5 . At the comoving scale l c ~ 10 kpc, Bg c ~ 10~ 50 G. In red (dots and solid) the spectrum 
Bg oc £~ 3 / 2 + u from inflation in our mechanism is shown, with and without preheating amplification. 
For this spectrum, Bg c ~ 10~ 29 G and Bi c ~ 10~ 34 G, respectively. We also show (blue dash- 
dots) the spectrum enhanced by helical turbulence (at l c ~ 10 kpc an enhancement of about 20 is 
obtained). This is to be compared with the dynamo bounds rescaled by a factor 5 x 10 3 (see main 
text) B seed > 2 x 10~ 27 G for a universe with critical matter density, and -B S eed ~ 2 x 10 34 G for 
a flat, low-density universe. 
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In Fig. |3| we display the spectrum Bi of the magnetic field for the different histories of 
amplification discussed in this section and compare with the comoving t c = 10 kpc seed-field 
bounds, which are -B see d ~2x 10~ 27 G in a universe with critical matter density (presently 
disfavoured by observations) and -B SC cd ~ 2 x 10~ 34 G in a flat, low-density universe domi- 
nated by dark energy The corresponding field strengths and bounds in a newborn galaxy 
can be obtained by multiplying with a factor 5 x 10 3 . 

To conclude, in this Letter we have presented a generic mechanism for production of 
gauge fields during inflation which predicts an almost scale-invariant spectrum for gauge 
fields, comparable with the scale-invariant spectrum of cosmological perturbations. When 
applied to electromagnetism, the resulting magnetic-field spectrum Bi oc r- 3 / 2 ~ r 1 can 
provide strong enough seed fields to explain the origin of galactic magnetic fields. 
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